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Abstract 



Using the q-version of the Darboux transform we obtain the general solution of g- 
difference Riccati equation from a special one by the action of one-parameter group. This 
. allows us to construct the solutions for the large class of g-difference Riccati equations 

' as well as g-diference Schrodinger equations. 



■ Introduction 

X 

■ In this paper we investigate the Darboux factorization method for g-difference version of 
Riccati and Schrodinger equations. It appears that this method which is by all means effective 
for differential Riccati and Schrodinger equation (@, ^ 3) leads to non-trival and interesting 
results in the g-deformed case too. Some of the new formulae have their undeformed version. 
They tend in the limit of g ^ 1 to the ones, which are well known in differential case. 

The Darboux factorization 
gives the well known correspondence between one-dimensional Schrodinger equation 



and the Riccati equation 



^^2'^{x) + V{x)]i;{x)=0 (2) 



-^u{x) = -u^{x) + V{x), (3) 
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where ^ 

u{x) = —lmjj{x). (4) 

This correspondence is a starting point for a search of the exact solutions of both equations 
above (see |^). 

We will investigate, in this paper, an analogue of the Darboux method for the pair of 
g-difference equation 

( i;{x) \_( R{x) S{x) \ ( \ 

and 

V{x) = dqu{x) — T{x)u{x) + R{x)u{qx) + S{x)u{x)u{qx). (6) 
whose solutions are related by 

It is clear that (|5|) and (J^) generalize (^) and (|3|) respectively. The Schrodinger and Riccati 
equations are obtained in the limit q ^ 1 under additional assumption that R{x) = = T(x). 
Let us recall here, that the g-derivative and g-integral are defined by 

d,^{x) = , (8) 

/ i;{x)d,t = Y,{l-q)q^x^{q^x) (9) 

•^0 n=0 

respectively, where < (7 < 1 . The standard derivative and integral are obtained for q = 1. 
However, the reason for the investigation of the (^-difference equations (|5|) and @ is not only 
that they generalize in a natural way the Schrodinger and Riccati equations. 

If one additionally assumes l — {l — q)xT{x) = in the real case and one takes q^ instead of 
the real argument x, the equation (|5|) appear to reduce to the three terms requrence equation 

V^„+2 = [1 - (1 - V-n+l + (1 - q?q^''^^S{q^+^)V{q^)^n, (10) 

for the function '■= "0(9") of the natural argument n £ IN U {0}. Hence, the g-difference 
equation (^) can be applied to those physical problems which are related to the theory of 
orthogonal polynomials (see |^) . 

The paper is organized in the following way. In Section 1 we will introduce the g-difference 
Darboux transform and we will integrate the equation (P) for the case of V{x) = 0. The 
action of g-difference Darboux transform will be presented in Section 2. There we find an 
one-parameter auto-Backlund transform for the g-difference Riccati equation and we show 
that it generates the general solution of (|6|) from a special one. Section 3 is devoted to the 
presentation of some extendend classes of solutions of the g-difference Schrodinger and Riccati 
equation. All results presented in Section 1 and 2 have well known differential counterparts 
and this aspect is also exhibited in the paper. 

1 The g-difFerence Darboux transform 

In order to solve the g-difference equation (|5|) by the iterative method we will rewrite it in 
the following form 

^-Af.)f 1, (111 



(p{qx) J \ 



2 



where 



A(x) = l-{l-q)x 



R{x) S{x) 



12) 



V{x) T{x) 

Let us assume here that 5'(x), V{x) and T(x) are continuous functions of a real argu- 

ment. Hence, the matrix sequence 

A(g"-ix) K{qx)K{x)=:K{x-q)n (13) 

is point wise convergent 

A(x;g)oo (14) 

to a matrix function A(x; g)oo. The inverse matrix function A.{x;q)^ is exactly the resolvent 
of equation (|5|), i. e. 

(:g)-<-)^'(:!S)' 

The problem of solving (^) thus is equivalent to the calculation of the infinite matrix product 



A(x;g)oo ■.= \{K{q''x). 



k=Q 



The above suggests the following transform: 

A(x) — > D{qxy^A{x)D{x) = A\x) 



(16) 



(17) 



^(x) 
fix) 



D{x) 



-1 



ip{x) 
ip{x) 



ip (x) 
(p'{x) 



where I?(x) is GL{2, M) - valued function of the real argument. It is obvious that the transform 



(|l^|l8|) preserves the form of the equation (|ll[) and the transformed resolvent A (x; 
related to the initial one by 

A'ix;q)^' = D{x)-'A{x;q)^'D{0). 



IS 



(19) 



Thus the virtue of the above transform is to find such matrix valued function D{x) for equation 
d^, which reduces unknown resolvent A{x;q)^ to some known one A (x;g)^^. 

We will find later the explicit form of the resolwent A{x;q)^ in the case when V{x) = 0. 
So, in order to integrate (|5|) it is enough to transform (|16|) to the uppertriangular matrix 
function A(x) by the use of (17-p^). Any matrix can be decomposed generically into the 
product of the uppertriangular and lowertriangular matrices. Thus, without loss of generality 
we can assume that 

-c(x) ? ) ^ 



D(x) 



After substituting ( |20D into we find that A (x) will be an uppertriangular matrix if and 
only if the function c(x) satisfies the q-difference Riccati equation (^). 
The g-difference Schrodinger operator factorizes into the form: 



+ V{x) = (dg + u{qx)) i-dg + ^.(x)) , 



(21) 



iff the function u{x) satisfies the equation 



is natural to call the matrix transform (|17 



6|) with R{x) = T{x) = and S{x) = 1. Hence, it 
I) the g difTerence Darboux transform. 
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We will use the identity 



n (1 - (1 - q)q''xf{q^x)) = exp ( [ 
n=o K^-lJo 

which is an easy consequence of the definition of g-integral. 



1 rinil-{l-q)tf{t)) 



dgt , (22) 



Proposition 1 IfR{x), S{x) andT{x) are continuous functions andV{x) =0 then 

B{x) 



A{x;q)c 







PXn ^ ln(l-(l-g)tr(f)) , 
'^-^f 1-g Jo t 



(23) 



where 



B{x) 



5(t) 



exp 



1 - (1 - g)ti?(t) VI 



1 /-^l^ l-{l-q)sR{s) 



.^" i-(i-g).r(.) '^^^)^^^ 



(24) 



Proo/; 

Since F(x) = we can assume that A(x; 5)00 is the uppertriangular matrix of the form 



n:Lo(i-(i-9)'?"^«('?"^)) B{x) 

n"=o (1 - (1 - g)'7"a;r(g"x)) 



where from the equation 

Hx; q)oo = A(gx; q)ooHx), 
we find that the function B(x) does satisfy 



(25) 
(26) 



B{x) = [!-(!- q)xT{x)] B{qx) - J|(l - (1 - q)q''+^xR{q'^+^x)) (1 - q) xS{x). (27) 

n=0 

The equation ( ^ ) is solved by iterative method. One gets finally: 

00 

B{x) = - JJ (1 - (1 - q)q''xT{q''x)) x 



n=0 



^1-(1 



(1 - q) q-^xSiq^'x) (1 - (1 - q )q^xR{q^x)) 

)q^xT{q^x)) 



„=0^ {^-q)Q''xR{q^x, 



(28) 



Substituting (E8) into (Eq) and using identity (p2D we obtain the formulae (p3D and (24) 



QED 



2 The solution of g-difFerence Schrodinger equation and auto- 
Baclund transform for g-difference Riccati equation 

We have obtained in Section 1, the resolvent function K{x]q)^ for the case of V{x) = 0, 
(see Proposition 1). Let us assume A (x) in ([l7|) to be uppertriangular matrix function. 
Applying the g-difference Darboux transform to A(x), with D{x) given by (pO[), we find the 
general solution of the g-differential equation (|5|) . 
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Proposition 2 Let us assume that functions R{x), S{x), T{x) andV{x) from the equations 
are continuous. Let ( iIjq{x) ipo{x) be a particular solution of one and put 

lPq{x) 



Uq{X} 



iI)q(x) 



Then the potential V{x) is given by 

V{x) = dqUo{x) - T{x)uo{x) + R{x)uo{qx) + S{x)uo{x)uQ{qx), 
and the general solutions of ^ is given by: 



il){x) = exp 



D + F 



\n{l-{l-q)t[R{t)+UQ{t)S{t)]} 



Sit) 



dqt I X 



1 



X exp 
(p{x) = F exp 
+uo{x) exp 



l-il-q)t[R{t)+uo{t)S{t)] 
l-q)s[R{s)+Uo{s)Sis)] 



* 1 1 
-In 



dnS d„t 



'^-qjQS l-{l-q)s[T{s)-Uo{qs)S{s)] 

'\n{l- {I -q)t[T{t)-Uo{qt)Sm 



l-q 



dqt I + 



1 f-\u{l-{l-q)t[R{t)+Uo{t)S{t)]} 



1-g 
D + F 



dqt I X 



S{t) 



l-{l-q)t[R{t)+Uo{t)S{t)] 



X exp 



*1 l.^l_q)s[R{s)+Uo{s)S{s)] 
l-qJo s'''l-{l-q)s [Tis) - Uo{qs)S{s)] ' 



The constants D and Fare related with the initial conditions by: 

D = ^{Q) , 



(29) 



(30) 



(31) 



(32) 



(33) 



Proof: 

In order to prove the formulae (^Tj- ^2|) and (|^) we assume in ( |l2| ) that A (x) is uppertriangular 
and we apply the g-difference Darboux transform ( p7| ) with 

D{x) = ( \ , ? V (34) 

The formulae (|30|) follows now from ([T7|) and from the property of A' (x) of being uppertri- 
angular. The formulae ( ^T| ) and (|32| ) are obtained from Proposition 1 and (|33|). 

QED 

As a corollary of Proposition 2 we obtained the auto-Backlund transform for g-difference 
Riccati equation. 
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Proposition 3 Let uq{x) be some special solution of the equation i^. Then the general 
solution of is given by 

u\x) = {B+uo) (x) = uo{x) + 



^ '^^P [T^ Jo y ^'^ l-{l-g)y[T(y)-uo{qy)S{y)\ "^iV ^ 

^ r gfa) ( I ry 1 l-(l-g)«[fl(s)+t.o(s)5(s)] , ^.J„/ ^ ' 

where t G iR. 
Proo/; 

The formulae ( |35| ) is obtained by substituting ( [3l| ) and ( ]3^ into (|^) and puting i = ^. 

QED 

The g-difference Schrodinger equation 

(-52 + V{x)) V(x) = 0, (36) 

is a special case of @ and is obtained by putting R{x) = T{x) = and S{x) = 1. From the 
Proposition 2 we may draw: 

Corollary 1 The solution of q-difference Schrodinger equation with the potential 

V{x) = dqUoix) + uo{x)uo{qx), (37) 

is given by 



]_ r ln(l-(l- g)tno(t)) 

i-qJo t 



il){x) = exp ( —- / dqt ) X 



'' + ''70 l-(l-.)Mt) I W io ^ l + (l-,)..o(..)'^V 

In the limit of g ^ 1 the g-difference equations (|5|) and (|6|) tend to their differential 
counterparts 

d ( i){x) \ _ ( R{x) S{x) \ ( iPix) 



dx V J V ^i^) T{x) J V ^(x) J ^^^^ 

and 

V(x) = —u{x) + {R{x) - T{x)) u{x) + Six)u'^{x), (40) 
dx 

where u{x) is given by (^). The Propositions 2 and 3 are valid in the limit g — > 1 too. 
Therefore, we can apply them to the differential case and reproduce in such a way the formulae 
for general solutions of (39) and (HO). They are given by 



i){x) = exp {R{t) + uo{t)S{t)) dtj x 

L> + F^ S{t)exp(^J^ [T{s) - R{s) -2uo{s)S{s)]ds^ dt 
ifix) = Fexp (^^ (T(t) - no(t)5(t)) dt^ + 
+uo{x) exp ( / {R{t) + uo{t)S{t)) dt] x 



(41) 
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D + F^ S{t)exp(^J^ [T{s) - R{s) -2uo{s)S{s)]ds^ dt 



(42) 



and by 



^ ' '^l + t J,^ Siy) exp {jy [Tis) - - 2uo{s)Sis)] ds) dy ' ^^'^ 



Here ( ^0(2;) ^^0(2;) ) and uo{x) are some special solutions ( |39|) and (|40|). 

In order to describe the properties of the family of solutions u^{x) = (^B^uq^ (x), t £ M, 
given by (|35|) let us formulate: 



Proposition 4 

• The transforms /orm one-parameter group 

BtoBZ=B+^,^, (44) 
which acts transitively on the space of all solutions of the q- difference Riccati equations 

• The solutions u^^{x), u^^{x), u*^(x), n*''(x) satisfy the unharmonical superposition prin- 
ciple 

{U'^{X)-U'-^{X)) {U'-{X)-U'^{X)) ^ (t4 - tg) (tl - ts) 

(^t3(2;)-n*i(x))(tX*2(x)-^/*4(x)) (t3_fi)(f2-t4)' ^ ^ 

/or ti, t2,t3,t4 G iR. 
Proo/; 

• According to Proposition 3 any solution of (||) is given by BgUo for some s G IR. Since 
Bt^ (BtUo) is a solution of @ we have 

Bt,oBluo = Btuo. (46) 



One can find s € iR by evaluation of both sides of (^6| ) at x = 0. Thus, using (33) we obtain 

ti + (t2 + ^^o(0)) = s + 'Uo(0). (47) 
• The equality is obtained from (|35|) by direct calculation. QED 
The right hand side of (^5|) is invariant with respect to the real fractional transformation 

"*^ + ^ ^^G5L(2,iR) (48) 



* cti + d' \ c d 

of the parameters ii, ^2, ^3 and t4. Hence, the left hand side of (|45|) is 5L(2, iR)-invariant 
too. 
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3 Some integrable cases 

In the Proposition 3 we constructed the transform which generated the general solution 
of the g-Riccati equation (^) from a given special one. It was done by the action of one- 
parameter group of transformations {B^}, t G M. The formula (|35|) which defines the group 
{Bf}te]R action does depend on the potentials S{x) and T{x). It does not contain the 

potential V{x). This allows us to define some method of creation of new integrable systems 
from a system which one knows how to integrate. In order to do this let us intoduce some 
notation. 

By T we will denote the map 

Zu{x) := -u{x). (49) 

The operator which acts on the function u{x) on the right side of g-Riccati equation (^) will 
be denoted by i.e. 

TZ^u{x) := dqu{x) — T{x)u{x) + R{x)u{qx) + S{x)u{x)u{qx). (50) 

By TZ- we will denote the operator 

TZ-u{x) := —dqu{x) + T{x)u{x) — R{x)u{qx) + S{x)u{x)u{qx). (51) 

It is clear that 

n+oI = TZ- and 7^_oT = 7^+. (52) 
Hence, if u{x) is the solution of the equation 

n+u{x) = V{x), (53) 

then Zu{x) does satisfy 

Ti-olu{x) = V{x) (54) 

and vice versa. The one-parameter groups ( p5| ) for the equations ( |53[ ) and (54) will be denoted 
by {S^lfgjR, and {B^}teM respectively. They are related by 

=Io e+ o I. (55) 

After aplication of the transform 

:= J o e+ o Z o S+ o . . . o Z o S+ o J (56) 

to uo{x) which one assumes to be the solution of (^) with the potential Vo{x) we obtain the 
solution 

u{ti, ...,tn,x) := B;'^ j-^uo{x) (57) 
of ( ^3|) with the some new potential 

TZ+u(ti, ...,tn,x)= V{ti, ...,tn,x), (58) 

which is n-parameter deformation of the initial one. The same function satisfies also 

TZ-U(ti, ...,tn,x)= V{t2, . . ■,tn,x). (59) 



In such a way we obtain the family of g-difference Riccati equations (|58|) and (|59|) generated 
by uo{x) and by the transform (^). 
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It is worth to mention here the group hke property of B, 



tl...tr, 



^o...o — 

^tl...t„ ° ^Sl...Sm ~ ^tl...tn-ltn+SlS2-..Sm' 

In particular of case n = 1, applying {B^}teiR to (|37|) we obtain 

TZ+u{t,x) = V{t,x), 

TZ-u{t,x) =1Z^uq{x), 
where u{t,x) = B^uq{x). Simple calculation gives a solution 

u{t, x) = Uq{x) — 

tew [j^ Jo y^^^T^iT^qj^^I^'iiy) 



(60) 



(61) 
(62) 



for the potential 



2dn 



y(t,a;) = yo(a:)- 
/ Pxn r ^ In i+(i-g)wfa) J „. 



(63) 



(64) 



1 + t f ^ 1 exD i-^ i In ±Hi:i9l£!M£l(7 s U ^ 

^^'■JO l+(l-g)j/«o(y) ^Vl-gJO s l-{l-q)suQ(qs)"'1^ ) 

In the limit of g ^ 1 the equations ( [6l| ) and (62) correspond to the proper differential equations 
and (|63) tends to the solution 



ii(t, rc) = uq{x) — In ( 1 + t [ exp (2 / tio(^)d2; ) 

ox \ Jo \ Jo ) 

of (^) with potential given by 

y(t, x) = Voix) - 2^ ln(^l + t exp (^2 no(2)dz^ dy 
Combining (^) and (|^) we find the formula 

u{t, x) = l\{l-q)xl {V{t, x) - 7^_no(x)) - 



(65) 



(66) 



1 2 



il-q)x-{Vit,x)-n^Uo{x)) 



+ 2{V{t,x)+n-Uo{x)) 



(67) 



which expresses the solution u{t, x) by the t-deformed potential V{t, x) and by the initial 
solution uo{x). 



Example 

The function 

uo{x) = ax", 

where a > — 1, a, a G iR, is the solution of (|37|) with the potential 

Vo{x) = a^-^x"-^ + a^g^x^". 
I - CI 



(68) 
(69) 
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Let 

oo ^ 

exp^(x) = Y: (70) 



be generalized exponential function (see j^) with 



1 _ rpOl+l 



Using the formulae (|63|) and (|6J) we find that 



^ expj;(a3;"+i) 



Uft X) = ax" expa(-ag°x"+i) . . 

^ Jo cxp^(-ag«y"+i) "^i/ 



is a solution of ( p7[ ) with potential x) given by 

exp^(ag"+ix"+i) 



2^. cxpfl(-ag«a:°^+i) 



. , cxp„(ag"+iy'^+l) , \ A exp^ (ag^t^+^y'^+O , \ 

^"^''Jo expfl(-aqaya+i) "gi/^ I ^ ^ V'' Jo exp^(-ag2«+iya+i) "g^/ I 

If, for example, a = 0, f = a and g — > 1 then the potential (^) reduces to Rosen-Morse 
potential 

a exp (-2ax) - 6 + exp (2ax) 

l/(o, xj=a 2 — • 

(exp (—ax) + exp (ax)) 
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